Spatial Reflection and Associated String Order in Quantum Spin Chains 
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We investigate spatial reflection and associated nonlocal order in spin chain quantum systems. 
The proposed string order parameters, e.g., reflected via operations of the spatial reflection or 
combinations of it with spin reflection, are able to characterize a variety of physical systems and 
allow us to gain renewed insights to the statistical mechanism underlying phenomena such as the 
Haldane gap and quantum phase transitions. Besides revealing further the potential application of 
the generalized parity symmetry in numerical algorithm, we build an explicit scheme to determine 
the symmetry and the related string order for matrix product states so that one can construct ansatz 
models with presumed properties. 

PACS numbers: 75.10.Pq, 03.67.-a, 64.70.Tg, 02.20. Qs 



Order parameter has been one of the most impor- 
tant concepts in condensed matter physics since Lan- 
dau's finding of it in describing continuous phase tran- 
sitions with spontaneous symmetry breaking. The non- 
local string order, originally introduced by den Nijs and 
Rommelsei, was thought of a hidden antiferromagnetic 
Neel order and exploited to understand the mechanism 
of the Haldane gap 2 . Despite that Haldane ground states 
of integer-^* Heisenberg chains possess only short-range 
spin- spin correlations, they could have nonzero string or- 
der parameters. The latter was found to be a common 
feature of the Haldane phase resulted from breaking of 
certain topological symmetry* and the phenomenon has 
been investigated for various physical models^. A typi- 
cal example is the valence bond solid (VBS) mode£ and 
its higher-S 1 generalizations^, in which the string order 
parameter can be explicitly worked out. 

To determine the hidden order for a general quantum 
many-body system is a problem highly nontrivial. In 
a recent literature^, the presence of den Nijs-Rommelse 
string order is demonstrated rigorously in relation to the 
existence of local symmetries within the framework of 
matrix product states (MPSs). The revealed connection 
is somewhat universal in the sense that the MPS formula- 
tion indeed offers a general mathematical representation 
for quantum many-body states^ ' 10 ' 11 . As a consequence, 
the problem of determining the den Nijs-Rommelse string 
order for a quantum system is recast as that of finding out 
possible local symmetries in the system. At this stage, 
it is of interest to build a scheme to identify the local 
symmetry for general MPSs, since for every MPS with 
a finite representative dimension an ansatz model with 
local interactions can always be constructed-, such that 
the given MPS constitutes its ground state. 

In this paper we investigate the spatial reflection and 
propose novel string order to characterize spin lattice sys- 
tems. Differing from conventional string operators utiliz- 
ing local unitary transformations, the present string oper- 
ator acting on spins in between two boundaries is defined 
by the spatial reflection or combinations of it with local 
spin reflection. By virtue of the MPS representation, the 
proposed string order is shown to relate intimately to the 
symmetry with respect to parity or generalized parity 



transformations. The quantity is then applied to char- 
acterize a sort of models including the VBS states and a 
spin-i system with quantum phase transitions. Similar 
to the parity symmetry, the generalized parity symmetry 
could also be applied to reduce computational costs in nu- 
merical algorithm of the density matrix renormalization 
group (DMRG)i2ii2. In addition, we build an explicit 
protocol to determine from first principles the string or- 
der for MPSs with representative dimension D = 2 and 
present an example of modeling ansatz with presumed 
properties. 

The primitive form of the string order parameter em- 
ployed here is proposed as 



SO£ = lim (S?V {l) S° 



(1) 



where Sf (a = x,y,z) denotes the spin operator of the 
ith lattice site and the parity operator acts on I = j — i— 1 
consecutive sites between the boundary spins i and j as 
pW|s i+ i • • • s i+ i) = \s i+ i ■ ■ ■ s i+ i). This string operator 
and its generalized form of Eq. (JTUJ) can be exploited to 
characterize the hidden order for infinite lattice systems 
just as the original den Nijs-Rommelse string order. 

To be specific, let us resort to the MPS formulation for 
translationally invariant spin chains 



l*> = 



1 



VTrE 



N 
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where = 1, • • • , d specifies the lattice spin degrees of 
freedom and {A s } is a set of D x D matrices param- 
eterizing the state |^). The transfer matrix E con- 
tained in the normalization factor is given as E — 
using further the notations, say, 

E S" = T, s ,s'( s \ Sa \ sl )( AS )* ® As \ the quantity SO P of 
the state \*f?) is shown to be 

Tr[Es9E v{ i ) Es< f E N - 1 - 2 ] 
SO£ = lim lim ' 1 „„' -, (3) 



l^oa N- 



TrE N 



where 



Ep W = [{E T -) l ] T \ E T * = J2( A T ® ( AS f- (4) 
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Note that the state \^f) is identified to be par- 
ity symmetrioi^ iff there is (^\V^ N) \^f) = 
Tr(E T2 ) N /TrE N = ±1. The latter means that 
for an infinite spin chain the largest eigenvalue 
X' m of E T2 should have the same modulus with 
the A m of E. This fact together with the relation 
limi-oo -Epw = A' m (|4>^)($^|) T2 lead to that the 
quantity SOp could be nonzero iff the state |\&) is parity 
symmetric and the corresponding value is calculated as 

SO P = ±(9%\Es?Q*Z)(*Z\) T >E ST \9*), (5) 

where |v?m ) and |3>m ) denote respectively the left 
(right) eigenvectors of E and E T2 related to A m . 

Let us address a concrete example of the VBS ground 
state of the AfHeck-Kennedy-Lieb-Tasaki models repre- 
sented by {A s } = {a z , \/2a + ,V2<7-}. In this case, E and 
E T2 are Hermitian thus their dual bases of left and right 
eigenvectors recover simply the conventional orthonormal 
bases. Their largest eigenvalue \ m — 3 and the corre- 
sponding eigenvectors are \^ m ) = ^(|00) + 1 11)) and 

|$ m ) = 4=(|01) - |10)). The string order SO P is worked 

out to be | SOp | = yg, which is independent of the ori- 
entation a = x,y or z owing to the SU(2) symmetry of 
the model. Moreover, the extended model of the SO(5) 
matrix product stat o 7 ' 16 specified by 

A" = a z (g>a z , A ±x = y/2a g ®a±, A ±2 = V2a±(g>I (6) 

has also nonzero SOp. The corresponding matrices E 
and E T2 are Hermitian again and have the largest eigen- 
values A m = 5 with eigenvectors as 

= l">' = ^(|14>-|41> + |32)-|23)). (7) 

In terms of the S z quantum number of the s — 2 lattice 
spins, the string order is obtained as SOp = — i. 

For further application let us consider a spin-^ matrix 
product system specified by A = ( ° °), A 1 =(* »). This 

is an ansatz modeli 7 - described by a parent Hamiltonian 
with three-body interactions and it undergoes a quantum 
phase transition at the point g = 0. The system is parity 
symmetric and the matrices E and E T2 have the same 
eigenvalues A± = 1 ± g. By definition, the string order 
parameter SOp of the ground MPS is shown to take a dis- 
crete form with SOp = g as g < and SOp = ~ ( ffpfj 
as g > 0. On the other hand, the system possesses also a 
local Zi symmetry with respect to the spin flip |0) <-> 
The transverse den Nijs-Rommelse string order associ- 
ated with the local unitary U = exp (iirS x ), namely, 

SO^ = lim (S?U® l S x ), l=j-i-l, (8) 

|j-i|-»oo 

is obtained as \SO x D \ = g/(l + g) 2 or SOf, = for 
cases g > and g < 0, respectively. Thus the quantum 



phase transition of the system crossing the point g = is 
marked distinctly by the non-analytical behavior of both 
these string order parameters. 

As the quantity SO p reveals the hidden order for states 
with parity symmetry, we show below that the state with 
generalized symmetry with respect to combination of the 
parity and spin reflection indicates another type of non- 
local order. For a simple example consider an MPS 15 
represented by {A*} = {(J ^),(^ J)}. The state 

is parity absent but invariant under the combined opera- 
tion V^Uf N , where U P is a spin flip operator specified 
as Up = |1)(1| + |2)(3| + |3)(2|. This sort of generalized 
parity symmetry is identified for MPSs as there is 

d 

{A*)* = Y t U i >{XAiX- i ), i = l,-..,d (9) 

and it is readily verified that for the above example 
there exists X = °). Indeed, the kind of sym- 

metry was already known in fermion systems in perfor- 
mance of the DMRG procedure^. Owing to the fermion 
anti-commutation rule, the Jordan- Wigner transforma- 
tion will result in a parity asymmetric Hamiltonian which 
is, however, invariant under the combination of the parity 
and a simple spin reflection. 

Accordingly, one is led to consider a string operator 

SOg P = lim (S?V®U$ l Sf), (10) 
|i— »|— >oo 

where V^Up 1 describes the combined transformation on 
/ = j — i — 1 consecutive lattice sites in between i and j. 
For an MPS with generalized parity symmetry specified 
by Eq. ((HJ), this string quantity is shown to be 

cnQ _ (^E s? V-\\^{^ n \) T2 VE s? \^ n ) 

^0 G p — , (11) 

where V = I®X, \^^ R ) are eigenvectors of E related to 
A m , and the relation lim / _ >oc (£; / ) T2 = ^n(l* m ) (*ml) T2 
has been applied to obtain the equality. The string order 
of the MPS mentioned above Eq. ((9]) is then derived as 
SOgp = gj and SOg P = — SO GP = gj. 

Revelation of the existence of the generalized parity 
symmetry in spin lattice systems is important in many 
aspects. For instance, it can be exploited to reduce com- 
putational costs in numerical algorithm. In fact, it has 
long been realized in the DMRG algorithm— that for 
systems with parity symmetry the density matrix of en- 
vironment blocks could be achieved via a simple reflec- 
tion on that of the system blocks. The various existence 
of generalized parity symmetries can also be utilized to 
achieve a factor ~ 2 speedup and the density matrices 
of system and environment blocks are now connected via 
the generalized parity transformation. Moreover, in view 
that an MPS can be viewed as a ground state of an ansatz 
system with short-range interactions, it is of interest to 
construct the parent Hamiltonian of the modeling sys- 
tem with certain presumed parity symmetry. The latter 
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implies that the derived system might possess nonlocal 
order associated with the specified string operator. 

We now focus on the problem of determining the gen- 
eralized parity symmetry for MPSs so as to construct 
ansatz lattice models with prescribed properties. The 
question arises as: given an MPS l^), how to certify the 
relation © exists or not and how to find out such Up and 
X if they do exist? To this end, we employ an expres- 
sion equivalent with Eq. ([9]) but in terms of the transfer 
matrix 



e t = {x* ® x)e(x* ® xy 



(12) 



In comparison with Eq. ©, this equation focuses solely 
on the invertible matrix X and the task reduces to dis- 
tinguish a particular S, among all those satisfying E T = 
SES -1 , that could be decomposed into S — X* <g> X. 
The difficulty of the problem comes from that the similar 
transformation S connecting E T ~ E is not unique and 
the complexity even increases with representative dimen- 
sions. Intriguingly, we show below that for cases of D = 2 
the problem could be resolved by virtue of a particular 
realization of group homomorphism between the direct 
product group of unimodular linear transformations and 
the four-dimensional complex orthogonal groupie. 

Briefly, it was shown in Ref. [l8j that in terms of a sort 
of so-called pseudo-orthonormal bases, the representative 
matrix of the operator X\ ® X2, where X\^ denote ar- 
bitrary linear transformation of SL(2,C), is complex or- 
thogonal and constitutes a group element of SO(4,C). 
The result leads to an intuitive perception that the 
pseudo-orthonormal bases could be utilized to resolve the 
factorization problem for the mentioned similar transfor- 
mation S. To be specific, the following simplest pseudo- 
orthonormal bases (the so-called "magic bases" , see Ref. 
fl9l ) will be employed 

\ei) = ^(|0Q) + |ll)),| e2 ) = ^(|00)-|ll)), 

N) = ^(|01) + |10)),|e 4 ) = ^(|01)-|10)).(13) 

Suppose that the transfer matrix has full rank and is 
expressed as E = Yli=i Xi\^f")(^f \ according to its 
spectrum structure. Its transposition is then formed 
as E T = Ei=i'M*i'*)(*f*l and the series of similar 
transformations connecting them are given by S(ki) = 
Si=i where the parameters fcj 7^ are 

to be determined later. The distinguishing problem of 
Eq. (fT2"|) could be resolved according to the follow- 
ing protocol: i) Express the operator S(ki) in terms 
of bases (fT3|) and denote the representative matrix as 
D^(fei) = (e M |S , (/c i )|e t ,); ii) Presume that D S (D S ) T = 
[D S ) T D s = I and resolve the set of ki if they do exist; 
iii) Factorize the operator as S(ki) = X\ ® X 2 accord- 
ing to the correspondence of the group homomorphism 
SO(4, C) - SL(2, C) ® SL(2, C) and verify if there exists 
X2 = X{. Once X is figured out, the spin reflection Up 
can be worked out easily from Eq. ©. 



It is worthy to note that a slightly modified scheme 
is capable to determine the local symmetry hence the 
string order SOd for MPSs, wherein the condition reads 
as£ A i = Y,j U ij {XA^X- x ). The problem then becomes 
to find out a decomposable transformation S — X* ® 
X from all those satisfying E = SES- 1 . The S here 
is expressed as S(ki) — J2i &i|*f )(*fl by invoking the 
spectrum structure of E. The distinguishing problem can 
be resolved following the same protocol described above. 

Below we present an example to illustrate the scheme 
and construct the interaction model with presumed sym- 
metry. Consider a family of MPSs \^{g)} with 



A 



1 

g 



,A+ = 




9 



1 

1 



(14) 



where g is a real parameter. The corresponding E has the 
largest eigenvalue Ai = |(fi + \/4 + Q 2 ) with £1 = 1 + g 2 . 
By comparing the spectrum of E and E T2 , one finds that 
the state is parity absent except cases of g = 0, ±1. To 
reveal the possible symmetry and hidden string order, it 
is instructive to specify the four left eigenvectors of E as 

l*i,2) =7l,a|00) + |*L) = |01) ± |10), (15) 



■ g 2 ± V4 + fl 2 ) . In terms of bases ([T 



where 7^2 = s(l - 
the transformation S(ki) — 5Zi=i ki\Li)(Li\ is expressed 
as D s {h) =(fi % ), in which B 2 ={ 2 " 3 2, ) and 



2fc 4 / 



1 2 



fe t (l + 7l ) 2 iki(l-jf) 



(16) 



By substituting it into the condition D s (ki)[D s (ki)) T 
I, one derives 



, _ v-7172 

K\ — — 5 , 

7f - 7172 



U = 



\/-7i72 
72 - 7172 



k a = h = ^. (17) 



Thus one can infer that the specified operator S(ki) falls 
into SL(2, C)(8>SL(2, C) according to the described group 
Homomorphism. Simply, the decomposition S(ki) — 
X\ ® X2 of this case can be obtained by expressing S(ki) 
in the conventional bases {|00), |01), |10), |11)} as it has a 
diagonal form S(ki) = diag{^2 ; 1, 1, fi - 5 }. The yielded 
Xi t 2 and Up satisfying Eq. © turn out to be 



Xi,2 = diag{n»,fi-»} J U P (g) 



(18) 



where the Pauli operator a is defined in a two-state space 
{|+), |— )} and the reflection is taken in the Bloch space 
along n = (cos 6, 0, — sin 9) with 9 = arctan^. As a 
consequence, the string order parameters SO^p of ^(g)) 
are figured out in Fig. 1. 

The parent Hamiltonian with fc-local interactions al- 
ways exists for a given MPS as k > 21og d Z? 9 . It can be 
constructed by a sum of positive operators supported in 
the null space of a Hermitian matrix y^L4^, where A is a 
d k x D 2 matrix with A%' ih = (A h ■ ■■A lk ) af3 . For the 
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FIG. 1: The transverse string order parameters SOq P (solid 
line in the left), SOq P (dashed line in the left) and the lon- 
gitudinal SOqp (right) in the state \^(g)) specified by Eq. 
(fT4)> . It is shown that there is SOcp(-ff) = -SO y GP {g). 



state \^(g)) specified by Eq. ([Llj) one gets a Hamilto- 
nian with k = 2, i.e., H = J2« h{i, i + 1). As g = 1, the 
detailed form of h is expressed as 

-^fSf+i - ^(A^.+i + + ilf, (19) 

where A = (S x ) 2 - (S y ) 2 , B = {S z f - S z and K = 
(S?S? +1 ) 2 + (S+S^) 2 + (SrSt +1 f. The state |*(1)> 
turns to be a unique ground state of the model by 



constructions^ and it has dual symmetries under the 
parity and local spin reflection Up (I) [cf. Eq. ([IS])]. 
Consequently, the system at this point has non- vanishing 
string order via both passages of the parity and the lo- 
cal unitary, e.g., they are obtained as SOp = |, and 

SO% = ^p, respectively. 

Summing up, the symmetry related to spatial reflec- 
tion has been investigated for quantum spin chains and a 
set of renewed string order was proposed. The quantity is 
shown applicable to characterize a variety of physical sys- 
tems owing to the particular role of spatial reflection in 
spin lattice systems. Furthermore, revelation of the gen- 
eralized parity symmetry is important in many aspects. 
Besides that it can be exploited to reduce computational 
costs in numerical algorithm, we build an explicit scheme 
to determine it for MPSs so that one can construct ansatz 
models with presumed properties. In addition, the string 
correlation was shown related to the concept of localiz- 
able entanglement^, a quantity reflecting the quality of a 
quantum channel. It would be of interest to explore fur- 
ther the role of the string order variants in the context 
of quantum information theory. 

The author thanks Peng Li for helpful discussions. 
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